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ABSTRACT 

We present a detailed observational analysis of how merger-driven turbulence may regulate the star-formation efficiency during galaxy 
interactions and set the initial conditions for the formation of super star clusters. Using VLT/SINFONI, we obtained near-infrared 
imaging spectroscopy of a small region in the Antennae overlap region, coincident with the supergiant molecular cloud 2 (SGMC 2). 
We find extended H2 line emission across much of the 600 pc field-of-view, traced at sub-arcsecond spatial resolution. The data 
also reveal a compact H2 source with broad lines and a dynamical mass M^y„ ~ 10' Mq, which has no observable Bry or A'-band 
continuum emission, and no obvious counterpart in the 6 cm radio continuum. Line ratios indicate that the H2 emission of both sources 
is powered by shocks, making these lines a quantitative tracer of the dissipation of turbulent kinetic energy. The turbulence appears to 
be driven by the large-scale gas dynamics, and not by feedback from star formation. We propose a scenario where the H2 emission is 
related to the formation of bound clouds through accretion. The kinetic energy of the accreted gas drives the turbulence and powers 
the H2 emission. Broad H2 line widths of-order 150 km s"', similar to the velocity gradient of the gas across the field of view, support 
this idea. Within this interpretation, the compact H2 source could be a massive cloud on its way to form a super star cluster within the 
next few Myr. This scenario can be further tested with ALMA observations. 

Key words. Galaxies: interactions, ISM, starburst, star formation 



1. Introduction 

Detailed observations of the nearby Antenn ae galaxy merger 
NGC4038/4039 at a distance of only 22 Mpc (Schweiz er et alj 
1200 8) provide important benchmarks for our understanding of 
how the large-scale dynamics of galaxy interactions trigger star 
formation on much smaller scales. Contrary to the 'canonical' 
scenario of merger-induced star formation , which predicts in- 
tens e starbursts in the nuclear regions (e.g.. lBarnes & HernquistI 
11996) . most stars in the Antennae are formed off-nucleus, in 
a heavily obscured, gas-rich, infrared-luminous region where 
both galaxies permeate each other, the "overlap region"(e.g., 
.Vigroux et al. 1996; Klaas et al. 2010). 

Most of the star formatio n in the overlap region of 
20 Mo yr ' jZhang et al.1 l2001h occurs m massive super-star 
clusters (SS Cs) with masses larger than 10"^ and up to a 
few 10*^ Mq (W hitmore & Schweizer 1995; Zhang & Fall' 1999^, 
IWhitmore et al. 2010). Detailed studies of SSCs constrain the re- 



cent star formation history of the overlap region. At typical ages 
of few lO*" yrs, many SSCs have already expelled most of the 
gas and dust from their immediat e surroundings, and they may 
be b ecoming unbound the mselves ('Fall et al. 2005 HMengel et alJ 
5005; Gilbert & Grahamll2007: .Fall et al. 201Q) . 

Our understanding of how the galaxy interaction may trigger 
gas collapse and the very early phases in the formation of SSCs, 
is much less devel oped and relies most ly on theoretical argu- 
ments. For example. lScoville et al.l (Il986l) suggested that the col- 
lision between the two galaxies could trigger collisio n s betw een 
pre-existing giant molecular clouds. i Jog & Ostriked (Il988l) ar- 



* Based on observations with the VLT/SINFONI, Program ID 383.B- 
0789, with the VLT/CRIRES, Program ID386.B-0942, and with the 
CFHT/WIRCAM, Program ID 09AF98 



gued that shock-heated ambient gas may enhance star forma- 
tion by inc reasing the pressure in embedded molecular clouds. 
iKeto et al.l 12005) discuss their high resolution interferometric 
observations of molecular gas and the formation of massive 
star clusters in the starburst galaxy M82 within this scenario. 
However, observations of the Antennae overlap region may not 
be accounted for by the collapse of pre-shock clouds. The large 
surface density of molecular gas and i ts fragmentation in com- 
plexes with masses of several 10^ Mq dWilson et al.ll2000l) . two 
orders of magnitude larger than masses of giant molecular clouds 
in spiral galaxies, are evidence for cooling and gravitational frag- 
mentation of the diffuse gas compressed in the galaxy collision. 

First attempts to quantify the impact of the interaction on 
gas cooling an d fragmentatio n have bee n made wi th numeri- 
cal simulations (ITevssier et al.ll20 10: Karl e Tani2011 l 2010), but 
their description of the energy dissipation of the post-shock, tur- 
bulent multiphase ISM is still schematic. The simulations fol- 
low the turbulent energy of the gas over a limited range of 
scales. The mechanical energy of the galaxy collision is ther- 
malized in large-scale shocks and radiated away by the shock 
compressed post-shock gas. Observations of the galaxy-wide 
shock in Stephan's Quintet dAppleton et al.ll2006t ICluver et"al] 
120 lOl) show that a dominant fraction of the collision energy is 
not thermalized in such shocks, but cascades to smaller scales 
into turbulent motion within molecular gas. This results from 
the clum py multi-phase stru cture of the pre-shock interstellar 
medium (iGuillard et al .12009?) . The fact that most of the collision 
energy is transferred to the molecular gas has consequences for 
the energy dissipation and the gravitational fragmentation of the 
post-shock gas, and thereby for induced star formation, which 
have yet to be understood. It is the motivation of this paper to 
investigate how gas compression and turbulence driven by the 
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large-scale dynamics of the two interacting galaxies trigger and 
regulate star formation in the Antennae overlap region. 

Bright H2 line emission powered by shocks (rather than star 
formation) has bee n identified in th e diffuse interstellar medium 
of the Milky Way ("Falgarone et al."2005i). as well as in a signif- 
icant number of extragalactic system s including many interact- 
ing ga laxies ( Appl eton et alj2006ciCluver et alj201^IZakamskal 
l2010h . This makes emission-line observations of warm H2 an 
interesting tra cer of the d issipation of turbulent energy in the 
molecular gas. lHaas et alj (l2005i) proposed, based on ISO obser- 
vations, that H2 line emission in the overlap region could probe 
large-scale shocks driv en into m olecular clouds. This was later 
called into question by Brandl et al. ( 2009 ), who found lower H2 
line fluxes with Spitzer-IRS, and suggested that the bulk of the 
warm molecular gas may be heated by star formation after all. 
However, at a spatial resolution of ~ 5", and a spectral resolution 
of R = 600 (~ 500 km s ') Spitzer-IRS does not aflow to study 
the molecular gas at the scales relevant for the fragmentation of 
the molecular gas, <100 pc and ~100 km s respectively. 

We use near-infrared spectroscopy obtained with the imag- 
ing spectrograph SINFONI and the echelle spectrograph 
CRJRES (both on the ESO Very Large Telescope) to quantify 
the dissipation of kinetic energy in the molecular gas on scales 
of few 10s to 100s of pc. To this end we observed the region 
around "knot 5" of iBrandl et all (12009^. the brightest H2 peak 
in the overlap region. This knot, which is near the massive 
and IR-bright cluster No. 28 of Whitmore et al. (2010), coin- 
cid es with the super-gi ant molecular cloud SGMC 2 observed 
by IWilson et al] OOcicih in CO. SGMC 2 is one of the most 
massive molecular clouds in the overlap region, with a mass of 

4 X 10^ Mq. Our pointing also coincides wi th knot K2a in the 
Herschel dust imaging of iKlaas et al.l (l2010h . Thus, overall, by 
selecting an area with particularly bright line emission of warm 
molecular hydrogen, we observed one of the main sites of molec- 
ular gas and star formation in the Antennae overlap region. 

The paper is organized as follows. The observations and data 
reduction are described in Sect. 2. In Sect. 3 we focus on the 
structure and kinematics of the warm H2 gas. We identify two 
components of H2 emission: diffuse emission extended over the 
600 pc wide field of view of SINFONI, and a compact source. 
We show that, for the compact and diffuse extended emission, 
the H2 emission is powered by shocks rather than UV heating, 
and estimate the bolometric H2 luminosity of each component 
combining Spitzer and SINFONI observations (Sect. 4). In Sect. 

5 we propose an interpretation of the observations which relates 
the H2 emission to the formation of bound clouds through gas 
accretion. In Sect. 6 we discuss the observations and our inter- 
pretation within a broader astrophysical context, highlighting the 
role of turbulence for the formation of super star clusters and 
the regulation of the star formation efficiency in galaxy mergers. 
Section 7 gathers the main conclusions. 

2. Observations 

2.1. VLT/SINFONI imaging spectroscopy 

Our analysis is based on observations carried out 
with the near-infrared integral-field spe ctrograph SINFONI 
(lEisenhaueretal.1120031: iBonnet etal.l 12004 ) on the ESO Very 
Large Telescope. Data were obtained in service mode under 
Program ID 383.B-0789 (PI Nesvadba). SINFONI is an image 
slicer with a field of view of 8"x8" and a pixel scale of 250 mas 
in the seeing-limited mode. The spectral resolution is R = 4000 
at A = 2.2 yum. We observed three regions in the Antennae 



galaxy merger in the ^T-band, one near each nucleus and one re- 
gion in the overlap region. These pointings were selected based 
on their bright pure-rotational H2 line emission detected with 
Spitzer/IRS. This paper presents an analysis of molecular emis- 
sion of the pointing in the overlap region. 

All data were taken in June and July 2009 under good and 
stable atmospheric conditions. We obtained a total of 1800 sec- 
onds of exposure time split into individual exposures of 600 sec- 
onds. We adopted a dither pattern where one sky frame was taken 
in-between two object frames to allow for an accurate subtrac- 
tion of the night sky. Data reduction was do ne with the s tandard 
IRAF tools to reduce longslit spectroscopy (lTodvllI993l) . which 
we extended wi th a set of SINFONI-specific IDL routines. For 
details see, e.g.. iNesvadba et alJ (120071 l2008h . The telluric cor- 
rection was done using bright (A" ~ 6 - 8 mag) stars observed 
at a similar air mass as our target. We used the light profiles of 
these stars to measure the PSF of our data, FWHM - ff.'l x 0'.'6. 
Comparison wi t h hig h-resolution HST/NICMOS imaging of 
IWhitmore et al.l (l2OI0l) shows that our PSF estimate is robust. 

2.2. CRIRES iiigii-resoiution spectroscopy and CFHT 
narrow-band imaging 

To complement our SINFONI observations we also obtained 
a high-resolution longslit spectrum with CRIRES on the VLT 
(program ID 386.B-0942; PI Herrera) and near-infrared imaging 
of the Antennae through the H2 1-0 S(l) and the Bry filters 
with WIRCAM at the Canada-France Hawaii Telescope (CFHT 
program ID 09AF98; PI Nesvadba). 

With CRJRES we obtained a total of 5 hrs of integration time 
of the H2 1-0 S(l) line at 2.12 //m in the rest-frame for one 
pointing centered on the compact molecular source (Sect. 13.31 ) 
and at position angle PA = 67°. Atmospheric conditions were 
good and stable. We nodded our target along the slit to ensure 
a robust sky subtraction. Data were dark-subtracted, flat-fielded, 
and sky subtracted. The wavelength calibration was done on tel- 
luric OH lines straddling the wavelength of the H2 1-0 S(l) line. 
Using a 0.4" wide slit, we reached a spectral resolution of 6 
km s"', a factor ~20 higher than the effective resolution of our 
SINFONI data. This allowed us to analyze the line profile of the 
compact molecular sourc e at a spectral resolution which is com- 
parable to that reached bv lWilson et alJ (12000.) for CO. 

We determined the astrometry of the near-IR CFHT images 
through cross-correlation with the 2MASS /T-band image of 
NGC4038/39. We give a short description of these data in Sect. 
13.11 They have also been used to improve the astrometric accu- 
racy of the SINFONI observations (Sect.lZ4ll. 

2.3. Ancillary data sets 

We complement our near-IR spectro -imaging data with 
CO(l-O) observations obtained with the Caltech Millimeter 
Array. C. Wilson generously shared her zeroth and first-moment 
maps of NGC403 8/4039 with u s. For a full di scussion of these 
observations, see IWilson et all (I2OO0I llool . Our SINFONI 
pointing in the overlap region coincides with their super-giant 
molecular cloud 2 (SGMC 2). Within the spatial resolution of the 
CO d ata, the peak of SGM C 2 coincides with the warm H2 peak 
#5 in IBrandl et all (llool . In our analysis we ad opt the mass, 
size and line width of SGMC 2 listed in Table 1 of IWilson et al.l 
(2000). 

In addition, we use p ublished resu l ts fr om Spitzer/IRS mid- 
infrared spectroscopy bv IBrandl et al.l (l2009i) . Our pointing cor- 
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Fig. 1. Left panel: Central regions of the Antennae galaxy merger including both nuclei. The K-band continuum, H2 1-0 S(l), and 
Bry line emission are shown in blue, r ed, and green, respectively. All data were taken with WIRCAM on the CFHT. Grey contours 
show the CO(l-O) line emission from I Wilson et al.l (2000). The white square r e presen ts our SINFONI pointing. The H2 1-0 S(l) 
emission has the same morphology as the Spitzer IR emission in iBrandl et al.l ( l2009l their Fig. 4). The H2 morphology does not 
follow the continuum morphology, however, the Bry does. Right panel: Zoom onto the CO map of the overlap region (top, the 
square shows our SINFONI field of view) and the line-free /T-band continuum morphology seen with SINFONI (bottom). North is 
up, and east to the left in all panels. 



responds to their Peak #5. In Sect. 14. 21 we will use the rotational 
H2 line fluxes from their Table 5 to compute the bolometric H2 
luminosity. The integrated S(2) flux measured with the Short- 
High (SH) module is a factor 1 .4 higher than that m easured with 
the Short-Low (SL) module. iBrandl et all (l2009l) attribute this 
difference to the different apertures used for the SH and SL flux 
measurements. Since the aperture of the SH module (4'.'4x 1 1'.'3) 
is best matched to the SINFONI field of view, we scale the S(3) 
flux, which was measured only with the SL module, by a factor 
1.4. 



3. Observational results 

In Sect. 13.11 we present the morphology and kinematics of 
the molecular and ionized gas. In the two following sections, we 
highlight the presence of extended H2 emission throughout the 
field (with broad FWHM- 150-200 km s '. Sect. and the 
discovery of a compact molecular source, with bright H2 line 
emission, which shows no Bry, nor continuum emission in the 
K-band (Sect.O- 



2.4. Astrometry of the SINFONI images 



The SINFONI field-of-view is too small to register these data 
directly with astrometric catalogues. Therefore, to ensure that 
the absolute positional accuracy of our data is better than 1", 
we collapsed the SINFONI cube over the wavelength range of 
our CFHT Ks-hand image, which we previously registered rela- 
tive to the USNO catalog of astrometric standards, and matched 
the coordinates of both images. This gives a positional accuracy 
of ~0'.'5 for the SINFONI data. The position of the super-star 
cluster, the brightest continuum source in the SINFONI field 
of view, is within ff.'S of archival HST/NICMOS images (taken 
through the F187N and F160W filters). Curiously, we find a 2" 
of fset relative to the co ordinates of the super star cluster given 
bv lWhitmore et al.1 (1201 a cluster #28 in their Table 8) 



3.1. Identification of individual components 

The left panel of Fig. [1] shows the overall morphology of 
the continuum-subtracted H2 1-0 S(l) and Bry line emission in 
the Antennae. H2 emission is found in the two nuclei and in the 
overlap region. Bry is found in the overlap region with a differ- 
ent morphology than the molecular gas. Our SINFONI observa- 
tions targete d the brightest peak in pure-rotational molecular line 
enussion bv lBrandl etaH(l2009h . The field-of-view of our obser- 
vations is marked by the white box. The right panel of Fig. [1] 
shows two zooms into this region. 

In the SINFONI data cube we identify several emission lines 
from molecular (H2 1-0 S(0)^S(3), H2 2-1 S(l)), and ionized 
(Bry, Hei) gas. Fig.|2]shows the maps of emission-line surface 
brightnesses, velocities, and FWHMs for H2 1-0 S(l) and Bry. 
A position-velocity diagram is shown in Fig. 1. Mean veloc- 
ity and line widths were measured with Gaussian fits, carefully 
masking bad pixels. We also constructed a line-free continuum 
map (lower right panel of Fig.[TJ by averaging the flux across the 
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Table 1. Parameters of the super star cluster in our field. 



Obs 


Size 


Age 


Ak 


Mass 


-^ion 


K 


Reference 




n 


Myr 


mag 


xlO' Mo 


xlO^- S-' 


mag 




HST 




4.8" 


0.8"* 


1.2" 






#28 from Table 8 in [1] 


Keck + VLT 


a'46" 


5n" 




4.1" 


1.3" 


15.2^ 


SSC-C in [2] 


VLT/SINFONI 


a'6« 


< 6 
4.5 


0.4 


6.0 


4.4"-'" 
6.6' 


15.0^ 


Our work [3] 



[1] IWhitmoreetaP ( l2010l) based on ACS and NICMOS observations, [2] [Gilbert & GrahamI j200D based on Keck/NIRSPEC-SCAM and 
VLT/ISAAC observations, [3] For details see Sect.lSTTI 

" from the comparison of integrated photometry UBVIHcy with population synthesis models, * estimated from measured E{B - V) = 2.20 and 
assuming Ry = 3.1 extinction curve ( Weingartner & Draine 2001), '" deconvolved FWHM from Lorentzian fit, based on the Bry equivalent width 
as an age indicator comparing it with Starburst99 (Leitherer et al. 2010) models, " not corrected for extinction, ' A'-band apparent magnitude for 
an aperture of 2", ^ It is unresolved. The angular size (FWHM) measured from a Gaussian fit corresponds to the size of the seeing disk, * from 
observed Bry luminosity, ' from the fluxes of the [Ne ii] and [Ne iii] lines at 12.81 and 15.56 yum . 
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Fig. 4. Integrated SINFONI spectra for the four different components in our pointing of the overlap region of the Antennae. All 
spectra are smoothed by 3 bins in wavelength. Around 2.0 //m the atmospheric transmission is very low so we did not plot this part 
of the spectra. We mark in red the wavelengths of the strong night sky lines seen in the spectra (OH lines). In the right panel, blue 
contours mark the aperture from which each spectrum was extracted. 



full band in each spatial pixel after masking emission lines, bad 
pixels, and wavelengths affected by strong night-sky lines. 

We identify 4 components in these maps: a compact contin- 
uum source (hereafter the super star cluster, SSC), two compact 
emission line sources (hereafter the molecular, and the ionized 
compact source, respectively) and diffuse extended emission. 



The three compact sources are labelled in Fig. |4] The SSC has 
previously been identified with HST and Keck. The main param- 
eters of the cluster are listed in Table [1] We estimated the stellar 
mass to 6 X 10'' from the A'-band and Bry flux es corrected for 
extinction (see Sect. 14.2.1b . using Starburst99 jLeitherer et al.l 
1999) . assuming solar metallicity, an instantaneous burst of age 
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Fig. 2. H2 1 -0 S( 1 ) and Bry morphologies and kinematics. North 
is up, and East to the right in all panels. Top: Line surface bright- 
nesses in units of xlO"'^ erg s"' cm"^. The contours show the 
K-band continuum in steps of 15%, starting at 10% of the peak 
intensity, 45 x 10"^" erg s"' cm"^ A"'. Center: Velocities in 
km s ' relative to the mean recession velocity in our field-of- 
view of 1556 km s Bottom: Measured FWHMs in km s 
All maps are spatially smoothed by averaging the original data 
over 3 pixels x 3 pixels (0'.'4 x 0'.'4 pix"'). The contours in 
the mid and lower images show the line surface brightness. 
Levels are 0.1, 0.2, 0.3, 0.5, 0.7, and 0.9 the peak intensity 
(3.8x10"'^ erg s"' cm^^ for H2 and 2.6x10"'^ erg s"' cm^^ for 
Bry). Offsets are relative to the peak in the K-hand continuum 
map, a(72000) : 12''01'"54!753, 5(72000) : -18°52'51'.'44. 



4.5 Myr, and a Salpeter initial mass function with lower and up- 
per mass cutoffs of 1 and 100 Mq, respectively. Our upper age 
limit of 6 Myr is based on the absence of CO band-heads in the 
K-band spectrum and is consistent with previous measurements 
(Table [T]). This cluster is one of the most massive clusters in 
the overlap region, and one of the brightest in the near-IR. The 
characteristics of the molecular and ionized compact sources are 
presented in Table [3] 

The integrated spectra of each of the four components are 
presented in Fig.|4] In Tabled we list line fluxes, mean veloc- 
ities and line widths for the extended emission, and the molec- 
ular and ionized compact sources. To construct the spectrum of 
the extended emission, we integrated over the full spatial ex- 
tent where H2 1-0 S(l) emission is observed, and subtracted the 
contribution from the compact sources. For the compact sources, 
we integrated within circular apertures of Cf.'G, (X'S, and Of.'l (for 
the stellar cluster, the molecular and ionized sources, respec- 
tively) matched to their sizes. For the molecular and ionized 
sources we subtract the nearby background using an outer annu- 
lus. For the SSC, we do not see a distinct emission-line compo- 
nent that could be separated from the surrounding line emission. 




1200 1400 



1600 



1800 2000 



Velocity [km/s] 



Fig. 3. H2 1-0 S(l) position - velocity diagram from north to 
south, at the position of the compact H2 source. The peak corre- 
sponds to the compact H2 source (Sect. [33). Positions are rela- 
tive to this source. 



Therefore, we did not subtract background line emission from 
the spectrum of the SSC. Had we done so, we would have ob- 
tained a pure continuum spectrum with no evidence of emission 
lines (not shown). 

The spectra suggest that the astrophysical nature of these 
sources must be very different. For the extended line emission, 
we observe emission lines from both molecular and ionized gas. 
The molecular source shows only H2 lines and the ionized source 
only the Bry and Hei emission lines. Fig. |5] presents a zoom 
onto the Bry and H2 1-0 S(l) lines for the extended emission- 
line region and the compact molecular source to illustrate the 
differences in line intensities and kinematics. For the extended 
emission, the velocities of the molecular and ionized gas are off- 
set by ~ 100 km s"' . Throughout the rest of the paper we focus 
on the molecular emission, i. e., the H2 extended emission and 
the compact H2 source. The analysis of the ionized gas will be 
presented in a future publication. 

3.2. Extended molecular emission 

Fig. |2] shows extended H2 line emission over most of the 
field-of-view, with a morphology that is very different from that 
of Bry which is more concentrated around the SSC. The velocity 
field shows a systematic gradient of ~ 200 km s"' from north- 
east to south-west over a projected distance of 600 pc, giving an 
average velocity gradient of 0.35 km s"' pc"'. Figure [3] shows 
the position - velocity diagram for the H2 1-0 S(l) line, along 
the direction of the velocity gradient (P.A. 15°). The diagram is 
centered on the compact H2 source and reveals the velocity gra- 
dient and the high turbulence in the field. The H2 line widths of 
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FWHM = 150 - 200 km s"' are comparable to the total gradient 
over the SINFONI field-of-view. The H2 line width is larger than 
those measured in Bry for the ionized gas (Fig.|5]l. 

T he H2 rotational lines observed with Spitzer jBrandl et alj 
l2009h are not spectrally resolved, the highest spectral re- 
solving power of these data is R - 600, coiTesponding 
to FWHM = 500 km s-i. The CO(l-O) velocity (Vlsr = 
1470 km s"') agrees with the mean, luminosity weighted, veloc- 
ity of the H2 extended emission (Table |2]i. CO shows a velocity 
gradient as large as that of H2 but over larger distances. CO ob- 
servations with a higher angular resolution are needed to resolve 
this gradient at the same scale as H2. The warm H2 gas traced 
by the near-IR emission lines has line widths more than twice as 
la rge as those meas ured in CO at the same position (73 km s ' 
in lWilson et ^1120001) . 




^ 



Extended Emission 



^.^Moleculor Source 



1000 1200 1400 1600 1800 
Velocity Observed [km/s] 



2000 



Fig. 5. Emission line profiles of molecular and ionized gas. Solid 
lines correspond to Bry and dot-dashed lines to H2 1-0 S(l). 
The kinematics of the extended molecular and ionized gas are 
different with an offset of ~ 100 km s The molecular source 
shows broad H2 line emission and no Bry emission. 



3.3. Compact molecular source 

The H2 flux map in Fig. |2] reveals a compact source with- 
out a counterpart in the ionized gas and continuum emission. 
The spectrum of this source in Fig. |4] and |5] is unlike that of 
the other emission components. To our knowledge, this is the 
first time that an extragalactic source with a A'-band spectrum 
showing only H2 lines is discovered. Similar spectra are seen in 
the Milky Way towards low-mass stars embedded in molecular 
clouds, which do not produce enough i onizing photons to po wer 
bright Hydrogen recombination lines dGiannini et al. I I2OO6I) . In 
these sources the H2 emission traces shock heating of H2 asso- 
ciated with the interaction of the stellar jet with the molecular 
cloud. 

Table |2] gives the H2 line fluxes and an upper limit on the 
Bry flux. The H2 1-0 S(l) flux of the molecular compact source 
accounts for 10% of the total emission from SGMC 2 in this 
line. To estimate the size of this source, we fit the azimuthally- 
averaged emission-line surface brightness profile along right as- 
cension and declination with Gaussian curves. The FWHM size, 
geometrically averaged over the two axes, is O'.'S. Correcting for 
the seeing we obtain an intrinsic source size of 0'.'5 correspond- 
ing to 50 pc. 




-400 -200 200 400 

offset velocity [km s"] 

Fig. 6. Integrated spectra of the compact H2 source observed 
with CRIRES. Fitting the data with a single Gaussian (green 
line) shows significant residuals. An adequate fit (red line) re- 
quires 2 Gaussian curves (blue lines), coiTesponding to a narrow 
and a broad component, respectively. 



We obtained high spectral resolution measurements of the 
line profile of the compact source with CRIRES (Fig.|6]l. Fitting 
this profile requires two Gaussian components with different 
line widths at a common velocity, with FWHM=50 km s ' and 
FWHM=160 km s"', respectively. We illustrate this by showing 
fits with one and two Gaussians in Fig. |6] In the two-Gaussian 
fit, the narrow component contributes 25% to the total H2 line 
flux. We verified that we recover the line profile measured with 
SINFONI after convolution with the line spread function of 
SINFONI. 

The compactness of the H2 source suggests that it is grav- 
itationally bound, and that we can estimate the gas mass using 
the virial theorem. For a homogeneous spherical cloud the bal- 
ance between gravitational and kinetic energy gives the virial 
mass as Mvir - 5 R cr^/G, where R is the radius, cr the ve- 
locity dispersion of the gas, and G the gravitational constant. 
In this formula, we estimate the diameter using the FWHM 
size of the near-IR H2 emission. We obtain a virial mass of 
1.3 X 10^ X (cr/21 km s"')^ Mq using the fine width of the nar- 
row component detected with CRIRES for cr. The CO data do 
not have the required angular resolution to separate the com- 
pact source from the extended emission. High angular resolu- 
tion observations with ALMA are required to estimate the virial 
mass from CO. The mass of the compact H2 source is a factor 
~ 30 smaller than the mass deri ved from CO obser vations for 
the SGMC 2 molecular complex (I Wilson et al.ll2000l) . From the 
virial mass and size, we derive the column density of the molecu- 
lar source and we obtain A^h = 6x10^^ H cm"^. This corresponds 
to a mean density of 6000 H cm"^. 

Radio continuum maps at 6 cm show a slight excess (50 juJy) 
at the position of our molecular source (see Figures 5 and 9 in 
lNeff&Ulvestadll2000h . Assuming that this radio flux is entirely 
thermal emission from ionized gas, it corresponds to a Bry flux 
of~ 5x10"'* ergs"' cm"^.This value is a factor 5 larger than the 
upper limit on the observed Bry flux in Table |2] after scaling by 
our estimate of the extinction coiTection (a factor 20, see Table[3] 
and Sect. l4.2.TI ). This slight discrepancy could indicate that the 
newly formed stars are concentrated at the center of the cloud. 
In this case their emission could be more attenuated than the H2 
emission. From the radio flux, we compute an ionizing photon 
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Table 2. Emission-line characteristics. 



Line 


Rest 
Wavelength 




Extended Emission 


Molecular Compact S' 


ource 


Ionized Compact S' 


ource 


Vlsr 


FluxxlO'^ 


FHWM° 


Vlsr 


FluxxlO'* 


FWHM 


Vlsr 


FluxxlO'^ 


FWHM 


H2 1-0 S(3) 


1.95756 


1490± 


:17 


6.1±0.4 


225±14 


1515+16 


6.5+0.2 


167+7 








H2 1-0 S (2) 


2.03376 


1461± 


18 


3.3±0.3 


234±18 


1538+17 


2.3+0.1 


140+11 








H2 1-0 S(l) 


2.12183 


1493± 


:17 


6.9±0.3 


215±9 


1547+16 


7.6+0.1 


146+2 


1481+18 


3.1+0.6 


<136 


H2 1-0 S(0) 


2.22329 


1465± 


:18 


2.8±0.2 


188±16 


1554+18 


1.6+0.2 


140+22 








H2 2-1 S(l) 


2.24772 


1450± 


18 


1.7+0.2 


146±20 


1530+26 


0.9+0.2 


243+44 




<0.36 




Bry 


2.16612 


1408± 


16 


9.5±0.2 


133±4 


1443 


<0.06 




1404+16 


25.8+1.1 


20+29 


He I 


2.05869 


1410± 


16 


5.2±0.2 


131+7 








1400+16 


16.3+1.0 


35+27 


Continuum* 








0.06 












0.1 





Wavelengths are given in yum, fluxes in erg s"' cm"-, and velocities and FWHM in km s"'. Error bars include the uncertainties in the wavelength 
calibration. 

" intrinsic FWHM after deconvolution with the instrumental resolution of ~136 km s"', * line-free /f-band continuum flux in units of 
erg s"' cm^^A"'. 



which, for a young (1-2 Myr) burst corre- 
sponds to a stellar mass of ~ 4 x 10^ M© based on Starburst99 
models. This stellar mass is 0.3% of the virial mass of the com- 
pact H2 source. 



Table 3. Observed characteristics of the compact sources. 



Table 4. H2 line emission: Spectral diagnostics 



Component 


H2 2-1/1-0 S(l) 


H2 1-0 S(l)/Br7 


Extended emission 






Full area 


0.2 


0.72 


Southern part 




3 


H2 compact source 


0.1 


>100 



Parameters H2 source H 11 source 



Size 


a; 5 


a'4 


e,. 


-20 




«e [cm^^^] 




58° 


Mass [Mq] 


1.3x10'' * 


4.8x10"*'' 


A'ion [phot, s"'] 


<2.6xl0'" 


8.5x10^° " 




2.6x10* ' 


5x10' « 



° assuming ionization balance, and not correcting for extinction, scales 

* virial mass estimated from the intrinsic 8(1) size and the CO line 
width. 

mass of the H 11 gas uncorrected for extinction, scales as e°-^ . 
from the upper limit of the Bry emission, uncorrected for extinction, 
derived from the Bry luminosity, scales as e^. 
^ H2 bolometric luminosity corrected for extinction (Sect. 3.3). 

* from the age and Nio„. Uncorrected for extinction, scales as e^. 



4. The nature of the compact and diffuse H2 
emission 

In this section we argue that, for both the extended emission 
and the compact molecular source, the H2 line emission is shock 
powered. We also estimate the bolometric H2 luminosity for both 
emission components. 

4.1. Spectral diagnostics: PDR or sliocked gas? 

In this section, we discuss the excitation of the H2 gas and 
the nature of the H2 emission for the extended emission and the 
compact source. We focus on the near-IR lines within the /T-band 
using spectral diagnostics that do not depend on extinction. 

The H2 excitation diagrams for the extended emission and 
the compact H2 source are presented in Fig. [8] We obtain two 



estimates of the gas temperature, the first by fitting the H2 popu- 
lation in the rotational states of the v=l vibrational state, and 
the second from the ratio between the J=l level of the v=l 
and v-2 states. For the compact H2 source we find 1700 K and 
2300 K, respectively, and for the extended emission 900 K and 
2700 K. The molecular gas has a range of temperatures. Most 
of the molecular gas is too cold to emit in the near-IR. The 
mass of warm gas seen in the near-IR is several orders of magni- 
tude smaller than th at observed in the mid-IR rotational lines by 
iBrandl et al.l (l2009l) . and also much smaller than the virial mass 
given in Sect. 13.31 However, even though the mass is small, the 
energies are very large as we will see in Sect. 14.21 

The ratio /?2-i/i-o = H2 2-1 S(l)/H2 1-0 S(l) is often 
used to discriminate between UV heating of H2 in photodisso- 
ciation regions (PDRs) and shock heating. The observed val- 
ues of the extended emission and compact molecular source 
are 0.2 and 0.1, respectively. W e compare the se values with the 
Meudon PDR code described in lLe Petit et all (|2006). Their Fig. 
27 shows that PDRs typically have ^2-1/1-0 = 0.5 - 0.6. Values 
of /?2-i/i-o in the range 0. 1 to 0.2 are only reached for very high 
densities and strong radiation fields (n > 10^ cm"^, x ^ 10^)- 
Such conditions do exist in massive star-forming regions, but 
only close to massive stars. This solution does not apply here 
because the compact H2 source is far away from the SSC, and 
has no strong ionized line emission indicating the absence of 
massive young stars. For the compact source PDR emission is 
als o ruled out by the ratio H2 1-0 S(0)/H2 1-0 S(l) (see Fig. 28 
of iLe Petit et~al]l2006h . The extended H2 emission-line region 
extends over the whole field-of-view and does not peak at the 
position of the SSC. We therefore discard UV heating and favor 
shocks as main heating mec hanism for the ex tended emission 
as well. Using the models of lKristensenI (l2007h . we find that J- 
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and C-shocks match the observed values over a wide range of 
densities (10'* - 10^ cm"^) and shock velocities (15 -50 km s"'). 

The ratio ^s(i)/Bi-y = H2 1-0 S(l)/Br7' is an additional, more 
empirical diagnostics to identify the origin of the H2 excitation. 
These diagnostics provide additional evidence of shock excita- 
tion of H2 in our data. We compare o ur results for the An tennae 
overlap region with data presented bv lP^ixlev et al.l (11990). Their 
data include 44 star-forming regions in 30 galaxies, where, typ- 
ically, the H2 1-0 S(l) line is weaker than the Bry line (only 4 
regions have a ratio higher than 1). Their flux ratios range be- 
tween 0.1 and 1.5 with a mean of 0.5 and a dispersion of 0.3. A 
few sources su ch as the mer ger NGC 6240 show a much higher 
7?s(i)/Bry value. [Puxlev et all (.1990) present models for a number 
of different scenarios and conclude that the line emission is most 
likely from Hii/PDR gas in massive star forming regions, for 
most of their sources. With our data, we obtain ^S(i)/Bry = 0.72 
and > 100 for the extended emission and the compact source, 
respectively. The compact source shows a very high ratio, which 
is far from the typical range in star-forming galaxies. The ra- 
tio of the extended emission is within the range of observed and 
modeled values by Puxlev et al. ( 1990). However, this is not true 
for all of the extended gas. In the region South of the bright Bry 
line emission (Fig.|7|l, we measure a ratio of 3. This value sug- 
gests that at least parts of the extended H2 emission is excited by 
shocks. 




Fig. 7. The color image shows the K-hand continuum, which 
peaks at the SSC. Pink and blue contours show the H2 1-0 S(l) 
and Bry emission, respectively. We also mark the compact ion- 
ized and molecular source (Hii source and H2 source, respec- 
tively). The green area corresponds to the gas which is likely not 
influenced by the SSC, we call it 'Southern part' . 



4.2. H2 bolometric luminosity 

Before discussing the H2 emission as a coolant of the molec- 
ular gas in Sect. |5] we need to estimate the H2 bolometric lu- 
minosity of the extended emission and the compact source. This 
will be used in Sect. l5.2l and Sect. l5.3l to discuss what may power 
the H2 emission. To achieve this we combine the Spitzer and 
SINFONI H2 line measurements. This is not straightforward to 
do because the near-IR H2 is attenuated by dust extinction and 
the Spitzer angular resolution is too low to separate the contri- 
butions of the extended and compact source to the mid-IR rota- 
tional line emission. First, we estimate the near-IR extinction for 



1 \ 

o Extended emission 
Compact source 




J I I I L 



7 8 9 10 11 12 
Energy [K] x 10* 

Fig. 8. H2 excitation diagram of the extended emission-line re- 
gion and the compact source (Sect. 3.1). The dashed and solid 
red lines show the gas excitation temperature derived from fitting 
the H2 population in the rotational states of the v=l vibrational 
state, and from the ratio between the J=l level of the v=l and 
v=2 states, respectively. 



both emission components (Sect. l4.2TTI ). Second, we estimate the 
bolometric H2 luminosities, under the assumption that the frac- 
tion of the H2 luminosity in the H2 1-0 S(l) line is the same for 
tiie extended emission and the compact source (Sect. l4rZ2b . 

4.2.1. Extinction correction 

To quantify the extinction in the extended gas, we compare 
the number of ionizing photons A^io,, obtained from the observed 
Bry emission line with that computed from the fluxes in the 
mid-I R [Neii]/ll2.8//m and [Neiii]/ll5.6//m fluxes dBrandl et al.l 
l2009h . The extinction corresponds to the ratio between the two 
estimated values of A^ion. This ratio provides the effective extinc- 
tion factor averaged over the entire field-of-view at the wave- 
length of Bry because the dust extinction of the mid-IR neon 
lines is much smaller than in the near-IR. The observed Aion 
in the field-of-view was measured by integrating the Bry flux 
over the entire SINFONI cube. Assuming case B recombina- 
tion, an electron density = 10^ cm"^ an d an electron tem- 
perature of = 10"^ K (see Table 8 of iHummer & StorevI 
1987), Aion = 4.4 X 10^^ phot. s"'. To co mpute the Mon from 
the Ne fluxes, we use the equations in lHo & Keto (2007) and the 
[Neii]/il2.8;um and [Neiii]/I15.6yum fluxes. We further assume 
that «e is smaller than the critical density,- 10^' cm"^, and adopt 
the galactic Neon abundance. We obtain Aion = 6.6 x 10^^ phot, 
s"'. Comparing with the result for Bry this suggests an extinc- 
tion factor of 1.5. 

The compact source has a very high column density. Ah = 
6 X 10^^ H cm"^ (Sect. [33T l. and therefore, the extinction must 
be much hig her. Using the Milky Way extinction curve for 
Rv = 3.1 in Weingartner & Draind (1200 ih this column density 
corresponds to Text - 30 in the A'-band. We roughly estimate 
the extinction correction from Text(^) with a simple toy model. 
We compute the attenuation due to the dust absorption for every 
point of a spherical and homogeneous cloud of constant density 
and line emission per unit volume. The flux that reaches the sur- 
face of the cloud depends on the path length I, and is attenuated 
by a factor e The optical depth t^(1) = K^p x I depends 
on the absorption coefficient at a given wavelength and the 
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mass density p. The value i s taken from the Ry - 3.1 curve in 
IWeingartner & Draind (12001 '). We sum the attenuated emission 
over the cloud and compare the result with the total emission for 
zero opacity. We find a correction factor - 0.7 x Text. We also 
run the model with a density and emission profile p oc r"^ for 
which we find - 0.6 x Text. In the following we only apply 
this extinction correction {e^ - 20) to the narrow component of 
the H2 line emission. In our interpretation of the data, the broad 
component is likely to come from the surface of the cloud, for 
which the extinction correction should be about a factor of 2, 
because we do not see the emission behind the cloud. 

4.2.2. Bolometric correction 

The angular resolution of Spitzer/IRS is only 5", too low 
to measure the total H2 luminosity, Hj"', including the pure- 
rotational mid-IR lines, of the compact source directly. To cir- 
cumvent this difficulty, we need to make an assumption. We as- 
sume that the H^"' /H2 1-0 S(l) ratio, after extinction correc- 
tion, is the same for the compact source and for the extended 
emission. The detailed shock modelling which would be needed 
to support or refine this assumption is beyond the scope of 
this observational paper The total H2 luminosity ( extended and 
compact source) is taken from the observations bv iBrandl et al.l 
( 1200 9). assuming that most of the H2 emission is from the first 
four rotational lines. Within these simplifying assumptions, we 
find bolometric H2 luminosities of Lh, = 5.4 x 10^ and 2.6 x 10^ 
Lq for the extended and compact emission, respectively. 

The extinction correction introduces significant uncertainties 
to these luminosities, in particular for the compact source. The 
total H2 luminosity (extended emission + compact source) does 
not depend on extinction because most of the H2 emission is 
from the mid-IR v=0 rotational lines measured by Spitzer for 
which extinction is negligible. It is the ratio between both lu- 
minosities which depends on the extinction correction. We set a 
lower limit of 8 x 10^ on the luminosity of the compact source 
by assuming that the extinction is the same for both components. 
The H2 1-0 S(l) emission represents 3% of the bolometric H2 
emission. 



5. From large scale gas compression to star 
formation 

The results presented in the previous sections raise three 
questions. How were the extended molecular gas component and 
the compact H2 source formed? Why is this gas so H2 luminous? 
What is the nature of the compact H2 source? In this section, we 
argue that the extended emission traces highly turbulent molec- 
ular gas formed where the galaxy interaction is driving a large 
scale convergent flow. We relate the H2 emission to the dissipa- 
tion of the gas turbulent kinetic energy and the formation of the 
SGMC 2 complex and the compact H2 source by gas accretion 
(Sect. 5.2). In Sect. 15. 31 we argue that the compact H2 source is 
a massive core on its way to forming a super star cluster. 

5.1. Gas compression and gravitational fragmentation 

The Antennae merger is close to peric enter passage when 
tidal forces are compressive (iRenaud et al.ll200 8, 2009). Thus, 
we consider that the velocity gradient observed with the H2 line 
emission across the SINFONI field is tracing a convergent flow 
driven by the galaxy interaction. The associated gas compres- 
sion has created the conditions for the gas to cool and to be- 



come molecular, hke in num erical simulations of convergent 
flows ("Hennebelle et al.ll200 8'). As observed in the galaxy col- 
lision in Stefan's Quintet, the mechanical energy of the interac- 
tion is not fully thermalized in large scale shocks dGuillard et al.l 
i2009 ). Due to the inhomogeneous, multi-phase, nature of the in- 
terstellar medium much of the gas kinetic energy decays from 
large to small scales before it is dissipated. Thus, we con- 
sider that it is the galaxy interaction that drives the molecular 
gas turbulence. This interpretation is supported by the fact that 
the large velocity gradient has the same magnitude as the tur- 
bulent line width, and by numerical simulations, which illus- 
trate how dynamical and thermal instabilities lead to the forma- 
tion of highly str uctured and turbulent molecular clouds where 
gas flows collide (jHeitsch et alj|2005l; IVazquez-Semadeni et al.l 
l2007llHennebeUe et al.ll200'8lK 

The gas surface density is high enough for the turbulent 
gas to gravitationally contract and fragment within the time 
scale over which th e tidal forces are compressive (~10 Myr, 
iRenaud et al.l l2008h . We compute the Jeans length (/Jjeans = 
5cr^/(27rGcr)) a nd mass ( Mjeans = ■^^i.,„.o') using the CO(l-O) 
observations of IWilson et al.l (12000) to estimate the gas veloc- 
ity dispersion (cr) and surface density (cr). The value of cr is 
30 km s '. The surface density cr is obtained taking the ratio 
between the total virial mass and the total area of the over- 
lap region cr ~ 500 Mq pc"^. We find 7?jeans - 350 pc and 
A/jeans - 2 X 10^ Mq. To compute the free-fall time scale Tfr 
= V37r/(32pG) we determine the mean density p from Rjenns 
and Mjeans- We find tj = 8 Myr. 

The Jeans length and mass are comparable to the size of 
the super giant molecular complex SGMC 2. It is thus likely 
that the complex has formed by gravitational contraction out 
of the gas compressed by the galaxy interaction. The gravita- 
tional contraction could contribute to the H2 velocity gradient. 
An estimate of the motions driven by the self gravity of the com- 
plex is obtained taking the ratio between the free-fall velocity 
Vff = 2 Rjeans/Tg - 90 km s"' and the Jeans radius. We find 
a value comparable to the velocity gradient measured with the 
H2 line emission (200 km s ' over 600 pc). This calculation 
provides an upper limit to the velocity gradient from the cloud 
contraction since it is unlikely that the complex is free-falling. 
The cloud contraction must occur at a slower rate than free-fall 
because it takes time to di ssipate turbulent energy dElmegreenl 
l2007llHufl^& Stahleiil2007l) . 

Thus, we propose that the SGMC 2 complex is formed by 
convergent gas flo ws drive n by the galaxy interaction and the gas 
self-gravity. Wils on et al.l (|2000) find that the gas mass inferred 
from the CO line luminosity matches the virial mass estimated 
from the the CO line width and conclude that the SGMC 2 com- 
plex is gravitationally bound. In the remainder of this paper we 
use 'CO complex' to refer to gravitationally bound gas in the 
SGMC 2 complex. The velocity gradient and line width of the 
H2 emission are both larger than the CO line width, suggesting 
that not all gas in the SGMC 2 complex is gravitationally bound. 
Most of the H2 emission arises from gas that is unbound. 

In the next two sub-sections we will interpret our H2 data 
within this scenario. We associate the H2 emission with the for- 
mation of the CO complex through gas accretion. We argue that 
the H2 emission traces the dissipation of kinetic energy, which 
is required for gas that is driven by the galaxy interaction to be- 
come gravitationally bound. 
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5.2. The energetics and formation of tine SGMC 2 compiex 

From observations of the Milky Way (iFalgarone et al .1120051). 
and, more recently, of extragalactic sourc es ( Guillard et al.l 
[2009; Nesvadba et al. 2010; Ogle et aT] |2010h . we know that H2 
line emission is a major coolant of the ISM associated with 
the dissipation of interstellar turbulence. The H2 line emis- 
sion may arise from shocks , as quantified by models such as 
iFlower & Pineau Pes Foreti (12010). and from friction betweens 
ions and neutral species in vortices (Godard et al. 2010). Models 
show that other emission lines contribute to the gas cooling but 
that they do not change the order of magnitude of the coo ling rate 
(iFlower & Pineau Pes For"itsll2010l: iGodard et al.ll2010h . Based 
on this earlier work, we consider the extended H2 emission from 
SGMC 2 as a quantitative tracer of the dissipation rate of the 
kinetic energy of the gas. 

Since the H2 line widths and the velocity gradient across 
SGMC 2 is twice the CO line width, the H2 emission does not 
come solely from the dissipation of the turbulent kinetic energy 
of the CO complex. A significant fraction of the emission must 
come from the dissipation of bulk and turbulent kinetic energy 
of gas driven by the galaxy interaction. Such a loss of kinetic 
energy is a required step for this gas to become gravitationally 
bound. 

In the following, we quantify this interpretation which con- 
nects the energetics of the H2 gas to the formation of the CO 
complex by gas accretion. We consider that the CO complex 
evolves in a quasi-static way assuming that virial equilibrium 
applies at all times. The virial equation relates the gravitational 
and turbulent kinetic energy, Eg^dv and Ziturb, and includes the 
external pressure Pext^ 



(1) 



where V is the volume. The turbulent energy is Eturb - 
3/2Mcocr^Q, where Mqo is the mass of the CO complex and 
(Tco is the velocity dispersion along the line of sight derived 
from the integrated CO spectrum of SGMC 2. For a fixed Pext, 
the exchange of energy associated with gas accretion and radi- 
ation is a ssociated with the de rivative of the enthalpy H of the 
complex dHuff & Stahledl2007h : 



V = -£tuA + 4 Pext V = £;„ - Lh, //h, (2) 



where is the energy input from gas accretion, Lh^ the H2 
luminosity, and /hj the fraction of the gas bolometric lumi- 
nosity that is radiated in the H2 lines used to compute Lh, ■ 
This fraction is smaller than 1 because some cooling occurs in 
lines which hav e not been measur ed. Combining H2 and far- 
IR observations. iMaret et al.l (l2009l) find that /hj is in the range 
~ 0.25 - 0.5 for shocks associated with gas outflow from low 
mass stars. We assume that the kinetic energy of the accreted 
gas is the main source of en ergy that balances radiative losses 
dKlessen & HennebelleHIoTO) . In doing this we neglect the en- 
ergy that comes from stellar feedback. The fact that the H2 emis- 
sion shows no enhancement around the SSC is an indication that 
stellar feedback does not have a significant contribution to the 

H2 emission (Sect. 16741 . 

We follow .Klessen & Hennebelld (l20Toh in introducing the 
efficiency factor e which represents the fraction of the accretion 
energy that drives turbulence. The turbulent energy dissipates 
over a time scale fciis(CO) - R/crco where R ^ 30 pc is the 
radius of the SGMC 2 complex (Mc Kee & Ostrikeiir2007l) . The 
balance equation between energy input and dissipation is: 



We compute the mass accretion rate necessary to drive the 
gas turbulence by equating the energy input and the radiative 
losses. This corresponds to a solution of equation (|2]i where the 
derivative of the gas kinetic energy /^turb and the term associated 
with V are both small with respect to the right-hand terms. Thus, 
we write: 



Ei„ = 3/2 Mace crj^ - LmJ/h. 



(4) 



(3) 



where cth, is the velocity dispersion along the line of sight de- 
rived from the H2 data. To quantify Mace and e, we consider 
only the turbulent component of the velocity field (FWHMh, ~ 
150 km s"'), assuming that accretion occurs from turbulent gas 
with a mean velocity equal to that of the CO complex. For cth^ = 
65kms"' and (Tco - 30 km s"', we find Mace = 5.3//h, M© yr"^ 
and e = 1.8 X /hj. For our interpretation to hold we must have 
e < 1 , and thus /h^ < 0.5. The fact that the H2 line width is 
larger than that of CO shows that this holds. In the following, for 
numerical calculations, we consider that /h, ~ 0.25. 

We scale Mace by the turbulence dissipation time scale to 
estimate the mass of gas from which accretion is occurring: 
Mdcc - M^cc X fdi,s(H2), where fdis(H2) - R/ctu^ is the turbu- 
lent dissipation time scale. The large scale flow can feed this 
mass reservoir, compensating for what is accreted by the CO 
complex, because the dissipation time scale is comparable to 
the dynamical time scale associated with the velocity gradient 
across SGMC 2 (3 x 10^ yr). We find Mace = 2.4 x 10V/h, M^. 
For /h, ~ 0.25, Maee is 25% of Mcq. Thus, the mean density 
in the CO complex may only be a few times larger than that of 
the accreted gas. This agrees with our estimate of the efficiency 
factor e. Based on their numerica l simulations and theoretical ar- 
guments, iKlessen & Hennebelld (2010) argue that the efficiency 
factor e is roughly equal to the ratio between the density of the 
accreting gas to the mean density of the bound system. 

Computing the ratio between the CO mass of SGMC 2 and 
the mass accretion rate, we obtain a time scale faee = 7 x 10^ x 

yr. For /h, ~ 0.25, this is slightly longer than the 10 Myr 
time scale over which the dynamics of the interaction may have 
been driving the convergent flow (Renaud et al. 2008). However, 
within the uncertainties in this rough calculation, we consider 
that the present accretion rate is close to the mean rate needed 
to account for the formation of the SGMC 2 complex as a result 
of gravitational fragmentation and gas accretion driven by the 
galaxy interaction. 

Our interpretation introduces a dynamical view of the 
present state of the SGMC 2 complex. The virialized CO com- 
plex is physically associated with gas which is too turbulent to 
be bound. This unbound gas is dynamically fed by the conver- 
gent flow. It contributes an accretion flow onto the CO complex, 
which has the necessary magnitude to drive the gas turbulence, 
i.e. to balance the dissip ation of the gas turbulent kinetic energy. 
iKlessen & Hennebelld (2010) have proposed a similar interpre- 
tation to account for the formation and subsequent growth of 
turbulent molecular clouds in the LMC. In Sect. 16.31 we discuss 
the impact of this interpretation on the star formation efficiency 
within SGMC 2. 



5.3. Tine nature of tine compact H2 source 

In this section we extend our interpretation of the extended 
H2 emission to the compact H2 source. Like for the SGMC 2 
complex we propose an interpretation where the H2 luminosity 
traces the dissipation of the turbulent energy of the gas and gas 
accretion. 
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As a first test whether this is a plausible interpretation, we 
compare the dissipation rate of the turbulent energy of the gas 
with the H2 luminosity of the source. Combining equations Q 
and we obtain an equation that we can use to estimate the 
mass of the compact source, Mdoud, from the H2 luminosity. 



3/2 Meioud xo-l/R^ex LnJfu. 



(5) 



where o-y = 20 km s ' based on the width of the narrow compo- 
nent in the CRIRES spectrum (FWHM = 50km s^'). We find 



= 2.7 X 10^ X e//H, X (Lh,/2.6 x lO'' Lq) Mq 



(6) 



This mass estimate depends on the ratio e/fn-.- Assuming 
that this ratio is ~ 1 as estimated for the CO complex, we find a 
mass slightly larger than the virial mass Mvir = 1.3 x 10^ Mq de- 
rived for the same velocity dispersion in Sect. 13.31 We conclude 
from this comparison that the H2 emission from the compact 
source may be accounted for by the dissipation of turbulent ki- 
netic energy if the source is a gravitationally bound cloud with 
a mass in the range 1 to a few 10^ Mq. With this interpretation, 
the compact H2 source is a gravitationally bound cloud, massive 
enough to form a super star cluster. 

Like what we have done in the previous section for the CO 
complex, we can use equation (HJi to estimate the mass accretion 
rate and thereby the cloud formation time scale. For ctu^, we use 
the width of the broad component (FWHM = 160 km s ') in 
the fine profile of the 1-0 H2 S(l) measured with CRIRES. We 
find: 



9.2 X 



0.25 



2.6 X lO^-Lo 



M0 yr 



(7) 



The ratio between Mdoud and Mace yields a cloud formation 
time scale by accretion of face = 2.9 x 10^ x (e/0.25)yr. This 
value is about twice the cloud dynamical time R/o-y. We would 
thus be observing a pre-cluster cloud early in its evolution. Even 
at this early stage it is remarkable that the fraction of the cloud 
mass that has been transformed into stars is very small (~ 0.3%, 
Sect. 13.31 ). We observe that the cloud has formed and become 
gravitationally bound without having formed massive stars. 

If our interpretation is correct, by observing the gas cooling 
through H2 lines, we may have discovered a massive cloud on 
its way to form a SSC within the next few Myr. However, this 
conclusion is only tentative, and will remain incomplete, until 
we obtain the missing information about the mass, the density 
structure and the kinematics of the bulk of the gas. The H2 line 
emission provides this information only for the fraction of warm, 
shock-heated gas. The missing information can be obtained with 
ALMA. 

6. Discussion 

Our observations are related to two main questions in the 
field of star formation. How do super star clusters form? What 
are the roles of turbulence and stellar feedback in regulating the 
star formation efficiency in galaxy mergers? We discuss how our 
observations and interpretation may be of general relevance for 
these questions. 

6.1. The formation of super-star clusters 

Star formation is the result of the hierarchical structure of the 
molecular interstellar medium estab lished by turbulent compres - 
sion and gravitational contraction (iMac Low & Klessenll2004|) . 



The available gas fragments over a range of masses and time 
scales. Massive clusters form at different times from the densest 
and most massive gas concentrations. For SGMC 2, this view is 
supported by the presence of two massive clusters (the near-IR 
SSC and the compact ionized source), which have formed before 
the compact H2 source. More clusters, too small to be identified 
individually, are likely to have formed or to be in the process 
of being formed. Thus, the formation of super star clusters in 
SGMC 2 can be understood within the same framework as the 
formation of smaller clusters in disk galaxies like the Milky Way. 

Our work highlights three physical parameters which may 
contribute to account for the formation of exceptionally mas- 
sive clusters in the overlap region of the Antennae. These may 
be of general relevance for the formation of super star clusters 
in other extragalactic objects: other interacting/merging galax- 
ies, and also starbursts and irregular dwarf galaxies ( Ket o et al] 
l2005l:IWeidner et al. 2010). (1) The presence of compressive mo- 
tion on large scaleswhich collects the gas. The formation of the 
SGMC 2 complex and that of the compact H2 source is driven 
by the galaxy interaction. The large scale dynamics triggers the 
gas compression necessary for their formation and their subse- 
quent growth by gas acc r etion. (2) As discus sed in earlier studies 
(Escala & Larsonll2008t IWeidneret aP 1201 0). shear is the sec- 
ond key parameter in the formation of massive clusters because 
it sets the maximum mass that the gas self-gravity can collect. 
In the Antennae, the present geometry of the inter action, which 
tempo rarily makes the tidal forces compressive (iRenaud et al.l 
l2008l) . is favorable to the formation of massive clusters. As dis- 
cussed bv lRenaud et"an (l2009l) such conditions occur repeatedly 
in galaxy mergers. (3) High turbulence is an additional key pa- 
rameter. Since the Jeans mass is proportional to cr* (see Sect. 
15.11 ). a high value of the velocity dispersion cr scales up the 
masses of the clouds formed by gravitational fragmentation. To 
form a super star cluster, it is not sufficient to bind a large mass 
of gas. It is also necessary that the star formation efficiency be- 
comes high when the gas mass is highly concentrated. It is likely 
that high turbulence is a key factor which prevents star formation 
to be efficient before the cloud mass has been concentrated by 
gravity. This is a prerequisite to form a den se, potentially bo und, 
cluster rather than a loose OB association ( lElmegreenl2008l) . 



6.2. The impact of turbuience 

The impact of turbulence on the star formation efficiency is 
a long standing topi c of research (Mac Low & Klessen 200#. 
iKrumholz & McKed (l2005i) and Pa doan & Nordlundl 020111) 
used numerical simulations to quantify the dependence of ffie 
star-formation rate per free-fall time on the Mach number of tur- 
bulence. Boffi agree that, independent of the mean cloud density, 
for a cloud near virial equilibrium, the fraction of ffie gas mass 
that is unstable to collapse is small. If this is correct, then the 
formation of dense clusters within a cloud supported by turbu- 
lence must occu r ove r several cloud c rossing ti mes as argued by 
iTan et alJ (l2006h and lKrumholz et alJ 02006) . Elmegreed (|2008[) 
expresses a diff'erent view point by proposing ffiat star formation 
becomes efficient within cores with a mean density > 10"* H 
cm"^ independent of turbulence. 

These ideas are being tested against observations of galactic 
star forming regions and cores, which are far more detailed than 
the present information we have on the compact H2 source. High 
angular resolution observations of molecular lines at mm wave- 
lengths, which unlike the near-IR H2 lines will trace the kine- 
matics of the bulk of the gas and the source density structure, 
would be needed for a detailed study of the formation of super 
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star clusters on the example of this source, and if it is consis- 
tent with present ideas based on observations of lower mass and 
less turbulent cores. Here we only make preliminary remarks 
that can soon be tested with ALMA. The mean density of the 
compact H2 source is ~ lO'* H cm"-'. We estimate its formation 
time scale to at least 10^ yr, but the star-formation efficiency is 
still very small (stellar to gas mass ratio ~ 0.3%). This remark- 
able result may indicate that the high turbulence (ID velocity 
dispersion ~ 20 km s"') has been very effective at preventing 
star formation. Star formation may occur rapidly once the rate 
of mass accretion becomes insufficient to drive the present am- 
plitude of turbulence. The fact that stars may be forming out of 
gas which has lost the turbulent energy it had during the initial 
gravitational contraction reduces the requirement on the star for- 
mation efficiency to form a bound stellar cluster. 

6.3. Star formation efficiency 

The star formation rate (SFR) in the Antennae merger has 
been estimated to 20 Mq yr"' for a total molecular gas mass 
~ 10'° Mq. Local values measured from CO and Ha images of 
the Antennae agree with the general coiTelation between the SFR 
and molecular gas surface densities, the Schmidt-Kennicutt law 
dKennicutt 1998: Zhang et al. 2001). The burst of star formation 
observed in the overlap region follows mainly from the high sur- 
face density of molecular gas with only a small (a factor ~ 3) 
enhancement of the star formation efficiency with respect to the 
Galactic value. This enhancement comes from the non-linear de- 
pendence of the SFR on the surface density of gas. The star for- 
mation rate per unit gas mass in the SINFONI field of view is 
similar to the global value. Based on their Spitzer and Herschel 
data, Brandl et al. (2009) and Klaas et al. (2010) find an SFR of 
~ 0.7 Mq yr"' for a molecular mass of4x 10^ Mq. 

To compute the star formation efficiency, we consider that 
star formation is triggered by gas compression and occurs over 
the characteristic ti me, ~ lOMyr, over which the tidal interac- 
tion IS compressive dRenaud et al.ll2008l) . This is also the cross- 
ing time across the SINFONI field of view for the CO line width. 
Over this time the gas mass converted into stars is 7x10^ Mq, 
which yields a star formation efficiency of 2%. The SSC ac- 
counts for a significant fraction of this stellar mass (see esti- 
mates of the cluster mass in Table[T]i. For the SGMC 2 complex, 
like for the compact H2 source, the low star-formation efficiency 
may result from the strong, driven turbulence. For the molecular 
complex, the star-formation efficiency may remain low until it is 
disrupted by the tidal interaction. If turbulence continues to be 
driven by on-going accretion, it is possible that the CO complex 
will be dispersed without losing its turbulent energy. It will be 
interesting to extend the present study to other complexes in the 
overlap region to test this idea. 

6.4. Stellar feedback 

It is interesting to compare this analysis and interpretation of 
observational data with numerical simulations. In their Figure 1, 
[Karl et al. (201 1) compare the star- formation rates in recent nu- 
merical simulations of the Antennae. All simulations predict a 
significant enhancement of the star-formation efficiency, relative 
to its value in the two spirals prior to interaction, after pericen- 
ter passage. In the simulations, the enhancement in the star for- 
mation efficiency is triggered by the gas compression and subse- 
quent gas cooling. For several runs the enhancement is one order 
of magnitude or more, i.e. larger than that derived from obser- 



vations. iKarTetal] dloTTl) present new simulations that quantify 
the impact of stellar feedback on the star formation efficiency, 
and argue that it is necessary to invoke stellar feedback to com- 
pensate for the gas cooling. Our H2 observations do not support 
this idea. 

In Sect. |5] we argue that the gas turbulence is driven by the 
galaxy interaction. Stellar feedback is energetically significant 
but not a key contributor to the turbulent energy of the molec- 
ular gas. This holds for the SGMC 2 complex and the com- 
pact H2 source. We start discussing feedback on the scale of the 
SGMC 2 complex using Starburst99 models. For a Salpeter IMF 
and continuous star formation over 10 Myr, the mechanical en- 
ergy associated with stellar winds and supernovae explosions is 
iMech(SF) ~ 4x lO^x (SFR/0.7 Mq yr"') Lq. Most of this energy 
is released by the stars in the S SC. The non -thermal component 
of the radio flux from the SSC (iNeff^ & Ulv estad 2000) indicates 
that the SSC has evolved to an age where the most massive stars 
are exploding as supernovae. For a stellar mass of a few 10^ Mq, 
the mechanical power from stellar winds and supernovae is also 
a few 10^ Lq. This value is one order of magnitude larger than 
the H2 luminosity of the extended emission, but the lack of en- 
hancement of the H2 line width towards the SSC is evidence that 
stellar feedback does not contribute significantly to driving the 
turbulent kinetic energy of the H2 gas. The energy from stellar 
feedback must be mostly transferred to the hot X-ray emitting 
plasma (Baldi et al. 2006). We plan to test this tentative conclu- 
sion with additional SINFONI data towards other super star clus- 
ters in the overlap region. 

For the compact source, we estimate the star formation rate 
by dividing the stellar mass of 4 X IO^Mq (from Sect. O by the 
formation time scale face ~ 3 x 10^ yr (from Sect. 15.31 1. We find 
SFR ~ lO^^Mgyr"'. For this rate, the mechanical energy from 
stellar winds is two orders of magnitude smaller than t he H2 lu- 
minosity of the compact source. Using formula (6) in iMatzneJ 
(120021) we verified that the contribution from proto-stellar winds 
around low mass stars is also much smaller than the H2 luminos- 
ity. 

7. Conclusions 

We presented an analysis of VLT/SINFONI near-IR imaging 
spectroscopy of the region with the brightest H2 rotational line 
emission in the Antennae overlap region, which has previously 
been identified with Spitzer/IRS spectroscopy. The region 
encompasses one of the supergiant molecular complexes, 
SGMC 2, discovered through CO interferometric observations. 
It corresponds to one of the brightest FIR knots in the over- 
lap region, and is near a young (~5 Myr) super star cluster 
previously identified in NIR imaging including HST imaging. 
Our imaging spectroscopy provides us with constraints on the 
spatial distribution, kinematics and excitation of the H2 and 
Hn gas out to a radius of 300 pc from this cluster. Based on 
these observations, we investigated how the large-scale gas 
dynamics may trigger the formation of the SSC and regulate 
the star formation efficiency. We have also discussed how the 
observations and our interpretation may be of general relevance 
for the formation of dense massive clusters and the efficiency of 
star formation in mergers. We list our main results. 

• We find extended near-IR H2 line emission across much 
of our SINFONI field-of-view with broad line widths 
(FWHM ~ 200 km s"'), larger than those measured from 
CO and Bry at the same position (which are 70 and 130 km s 
respectively). The line width is commensurate with a large-scale 
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velocity gradient across the field. We argue that this extended 
emission component traces a convergent turbulent flow driven 
by the galaxies interaction. Spectral diagnostics show that the 
H2 emission is shock powered and traces the dissipation of the 
gas turbulent kinetic energy. 

• The data reveal a compact H2 source (50 pc diameter) with a 
A'-band spectrum showing only H2 line emission. The H2 lines 
are spectrally resolved with a width ~ 150 km s"' (FWHM). 
The H2 emission from this source is also shock excited. The 
cloud virial mass is ~ 1 X 10^ Mq. The absence of Bry emission 
and of some obvious counterpart in the radio continuum set 
a low limit on the mass of newly formed stars (stellar mass 
fraction ~ 0.3%). To our knowledge, this is the first time that an 
extragalactic source with such characteristics is identified. 

• The width of the H2 spectra show that the SGMC 2 complex 
and the compact source are both associated with gas which is 
too turbulent to be bound. We argue that this suggests that the 
H2 emission is powered by gas accretion. We show that the 
required accretion rate is of the right order of magnitude to drive 
the gas turbulence and to account for the formation of both the 
SGMC 2 complex and the compact source by accretion. 

• By observing gas coohng through H2 fines, we may have 
discovered a massive cloud on its way to form an SSC within the 
next few Myr. However, this conclusion is only tentative, and 
will remain incomplete, until we obtain the missing information 
about the mass, the density structure and the kinematics of the 
bulk of the gas. The H2 line emission provides this information 
only for the warm shock excited fraction. The missing informa- 
tion can be obtained with ALMA. 

• We propose that the strong turbulence observed in the 
extended emission and the compact H2 source is of general 
relevance for the formation of SSCs in two ways. A high value 
of the gas velocity dispersion increases the masses of the clouds 
formed by gravitational fragmentation. It may also prevent star 
formation to be efficient before the cloud is fully formed, i.e. as 
long as turbulence is driven by accretion. Clusters may form out 
of gas which has lost much of the turbulent energy it had during 
gravitational formation of the pre-stellar cloud. If this is correct, 
then it reduces the requirement on the star-formation efiiciency 
to form a bound stellar cluster. 

• Our observations highlight the role of merger-driven turbu- 
lence in regulating the star-formation efficiency in interacting 
galaxies. Within the field of view of our SINFONI observations, 
about 2% of the molecular gas has been turned into stars. We re- 
late this small efficiency to the fact that the gas turbulence does 
not dissipate because it continues to be driven by on-going ac- 
cretion. 
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